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When a single two-level atom interacts with a pair of Laguerre-Gaussian 

beams with opposite helicity, this leads to an efficient exchange of angular 

momentum between the light field and the atom. When the radial motion 
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• is trapped by an additional potential, the wave function of a single localized 
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atom can be split into components that rotate in opposite direction. This 
suggests a novel scheme for atom interferometry without mirror pulses. Also 
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atoms in this configuration can be bound into a circular lattice. 
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I. INTRODUCTION 



It is well-known that light may carry both angular and linear momentum. When a light 
field interacts with matter, exchange of momentum and angular momentum between light 
and matter can occur. Laguerre-Gaussian (LG) light modes are known to carry orbital 
angular momentum. If one employs the paraxial approximation for the light field, simple 
expressions for the field amplitudes and its average angular momentum can be derived [1]. 
An easy way to produce such beams is using spiral phase plates [2]. 

Another important question is the separability of the total angular momentum into 
'orbital' and 'spin' parts [3]. The orbital part is associated with the phase distribution of 
the light field, and the spin part is connected with its polarization. This question is essential 
in the context of momentum transfer from light to the atom when one includes atomic 
internal degrees of freedom. It has been shown that 'spin' and 'orbital' angular momentum 
of the photon are transferred from the quantized light field to, respectively, the internal and 



1 



the external angular momentum of the atom. The interaction with a LG mode is a possible 
way to entangle internal and external degrees of freedom of an atom [4]. The transfer of the 
angular momentum of light to particles has been also experimentally demonstrated in [5], 
where trapped massive particles are set into rotation while interacting with the light field. 
Other authors have studied the cooling properties for atoms using LG beams [6]. Also LG 
beams have been proposed as a 2D trapping potential for Bose condensates [7]. 

Whereas angular momentum exchange between light and matter is a relatively new 
topic, the linear momentum exchange is a well-established issue [8]. It is well-known that 
two counterpropagating waves lead to a more efficient exchange of linear momentum between 
an atom and the light field than a single travelling wave. Using quantum language for a 
classical light field, one can describe such an interaction as a sequence of successive single 
photon absorption and emission events. This suggests that one may expect more efficient 
angular momentum exchange between a light field and an atom if one uses two LG modes 
with opposite helicity, e.g. counterrotating waves. 



We start with radiation fields propagating along the z-axis with wave number k and car- 
rying orbital angular momentum (Laguerre-Gaussian beams). If one considers the paraxial 
limit of these waves, the expressions for the light fields are particularly simple [1] 



where p, z, 4> are the cylindrical coordinates, u is the frequency and the integer index I is 
indicates the helicity of the LG beam. For two Laguerre-Gaussian beams with opposite 
helicity, namely I and —I, the total field can be written as 



We indicated already in the Introduction, that one expects a more efficient exchange of 
angular momentum between the light field and the atom in the configuration (2) than in 
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II. GENERAL FRAMEWORK 



E (p, z, 0, t) = E (p, z) <;'(**-<•*+'*) + c.C, 



(1) 




(2) 



a single LG mode. This expectation is based on the corresponding situation of momentum 
exchange between an atom and a standing light wave. In addition to the light field (2), 
the atomic motion in the radial direction is assumed to be confined by an extra trapping 
potential U (p) with cylindrical symmetry. 

The z— dependence of the amplitude E (p, z) is slow and can be ignored. Properly 
shaping the LG mode, the radial dependence of E (p, z) can be ignored on the characteristic 
width of the trapping potential U (p). Thus, we assume that E (p, z) ~ E is constant. For 
a two-level atom the Hamiltonian in the rotating-wave approximation can then be written 
as 

H = H + U(p) + 2tiuj R cos {e l(kz ^ l) \e) (g\ + e'^'^ \g) (e|) , (3) 

where ojr is the Rabi frequency of each of the travelling waves that create the standing wave, 
uj is the laser frequency, and 

Ho=2M + ^(\e)(e\-\g)(g\) (4) 

is the Hamiltonian for a free atom, with P the momentum operator of the atom, \g) and 
|e) indicate the ground and excited states, and u = [E e — E g )/h defines the transition 
frequency of a free atom. 

The dynamics of the atom is rather simple if the laser is far detuned. We assume that 

|A|»w fl , (5) 

where the detuning A is defined as A = u — uj. For an atom in the ground state, the 
excited state can be adiabatically eliminated, which leads to an effective Hamiltonian in the 
well-known form 

H=^ + U(p) + V(<j ) ), (6) 
where the light-shift potential is specified by 

V{<t>) = -ftficos 2 Z0 (7) 



III. TRAPPING IN COUNTERROTATING FIELDS 



Some general conclusions on the bound states of the Hamiltonian (6) directly follow from 
its symmetry properties. We introduce the unitary translation operator T defined as 

T\<P) = \<P + j), (8) 

where |0) indicates the states with fixed azimuthal angle. Since the Hamiltonian (6) is 
invariant for rotation about an angle n/l, it follows from a rotational version of the Bloch 
theorem that the eigenstates of this Hamiltonian are also eigenstates of T. The eigenvalue 
relation can be expressed as 

T |^. = e-^|^>., (9) 

where q is referred to as angular quasimomentum and j identifies the energy band. We 
consider a single energy band, and we suppress the index j. We can restrict q to the first 
Brillouin zone given as 

-l<q<l. (10) 

The eigenstates \^ q ) should be periodic in <p with period 2n, because a rotation over 2-n 
must leave the wave function invariant. The finite range of leads to a discretization of 
angular quasimomentum. On the other hand, a rotation over 2n is equivalent to the action 
of the operator T 21 . Since it follows from Eq. (9) that 

f 2l \q q ) = e- 2iwq \^ q ), 

we conclude that the only possible values of the angular quasimomentum are determined 
from the condition 

e -2Mrg = L 

Hence q must be integer, and each band contains 21 Bloch states. For example, for I = 
2 the first Brillouin zone contains only the four values q = —2,-1,0,1 of the angular 
quasimomentum. 



Also, in analogy to the case of an infinite linear lattice, one can introduce localized 

Wannier states |0 n ) in the usual manner, as Fourier transforms of the Bloch states 
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Obviously, the number of Wannier states within an energy band is equal to 21, just as the 
number of Bloch states. 

In Fig. 1 we plot the trapping potential (7) V (x,y) /Kl for I = 2,4 in Cartesian 
coordinates. When the potential is sufficiently deep, atoms can be bound in the angular 
wells, and the Wannier states are confined to a single well. An additional confining potential 
U (p) is required to trap particles in the radial direction, and to avoid their escape. Then the 
potential (7) can create a circular lattice, where particles are located near the minima of the 
periodic potential. A circular optical lattice has many applications, as discussed recently by 
several authors [11]. 

IV. DIFFRACTION IN COUNTERROTATING FIELDS 

Since the potentials have a cylindrical symmetry, it is convenient to express the kinetic 
energy in cylindrical coordinates, and we write 

P 2 h 2 ( d 2 Id d 1 d 2 \ 

2M ~ ~2M \lh 2 + ~ P ~d~p P lT P + ~fd^J ' ( ' 

The dynamics along the z axis is completely free. For simplicity, we assume that the radial 
potential U (p) is narrow, so that the radial motion is restricted to a ring with radius p , and 
we ignore radial dispersion in the present Section. We return to it in Sec. VI, where the 
effect of the radial dispersion is estimated. The motion of an atom in the angular direction 
is then described by the one-dimensional Hamiltonian 

H = --—-hncos 2 l<j ) , (12) 

which has the azimuthal angle as the only coordinate. The quantity / = Mp\ is the mo- 
ment of inertia. This Hamiltonian is the circular counterpart of the Hamiltonian for simple 



linear diffraction. The main difference is that the coordinate is periodic, which forces 
the angular wave number / to be integer. Diffraction of a single atom described by such a 
linear Hamiltonian has been extensively studied theoretically and experimentally by several 
groups [8]. 

Just as is usually done for linear diffraction, we consider the situation that an initially 
localized atom interacts with the optical potential during a small interaction interval [— r, 0] , 
where the atom picks up momentum from the lattice. The transition from the near field 
immediately after the interaction and the far field is described by free evolution. We assume 
the atom to be initially in its ground state and situated in a small segment of the ring. 
Since the angular wave function $(0) of the atom must be periodic at all times, we cannot 
represent a localized wave packet by a Gaussian. The initial state at the beginning of the 
interaction interval is taken as 



The state (13) can be written as a Fourier series, which is just an expansion in the angular- 
momentum eigenstates. This gives 




(13) 



with N to be a large natural number, and Cat is the normalization constant 
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(14) 




(15) 



with 



2N 



\ 



1 



(16) 




N + m 



J 
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The initial state (13) is localized around = 0, which is clear from the asymptotic form 



N 



cos 2 ^expj-^ j, (17) 

for large N. The half width in the azimuthal angle is of the order of y/2/N. From the 
asymptotic form of the binomial coefficient 



yN + m j 



, x :i / in' 1 



~ 2 ^ — = exp . . . 



we find the asymptotic expression of the Fourier coefficient 

2 V' 4 I m' 



exp -7T'- (18) 



This demonstrates that the half width in angular momentum is of the order of y/~N/2. 

If we take the duration r of the light pulse short and the moment of inertia / is large, so 
that h 2 l 2 r/ (21), no propagation occurs, and the kinetic-energy term can be neglected during 
the interaction. This is the equivalence of the standard Raman-Nath approximation applied 
by Cook et al [10]. Then the final state at time after the interaction is 

$ (0, 0) = $ (0, -r) exp (iQr cos 2 l<f>) . (19) 

This state can be expressed as an expansion in angular- momentum eigenstates, in the form 
Fourier series, which is just an expansion in the angular- momentum eigenstates. This gives 

$(0,O) = -^=]TC m e^, (20) 
y zt[ — 

m 

where 

C m = exp (int/2) i n 1Pm-2nlJn (^/2) , (21) 
n 

in terms of the ordinary Bessel functions. 

States with large angular momentum \m\ > N are initially not populated, whereas all 
angular momentum states get populated after the interaction. Thus, the configuration with 
two LG modes leads to more efficient exchange between the light field and the atom than 



a single LG beam. The physical interpretation is the same as for diffraction in the field of 
classical counterpropagating waves: an atom picks up a photon from the light beam with one 
helicity and emits a photon into the opposite one. In figure 2 we present a typical diffraction 
pattern calculated for the case that I > N. More precisely, we plot the angular-momentum 
coefficients \ip m \ 2 before the interaction, and the coefficients |Cm| 2 after the interaction with 
the circular lattice, for Qr = 6. In the latter case, the momentum peaks correspond to 
different values of n. The distance between neighboring peaks is equal to 21. The half width 
of each peak is of the order of y/N/2. 



As shown above, the angular-momentum distribution of an atom after the interaction 
with a pair of counterrotating LG beams can be broad. However, as a result of the Raman- 
Nath approximation, the angular distribution of the atom has not been modified during the 
interaction, so that \&(<f>, — r)| 2 = |<3>(0, 0)| 2 . In this chapter we investigate the spatial form 
of the atomic distribution in the far field, i. e. after free evolution of the atom over the 
ring. As before, the motion along the z axis is completely free and the radial motion is 
restricted on a ring. The initial state of this free evolution is given by Eq. (19), with the 
expansion in angular-momentum states given by Eqs. (20) and (21). For positive times, the 
atomic motion is still restricted to the ring with radius po by the confining potential U(p), 
and the evolution of the angular wave function $(0) is governed by the Hamiltonian (12) 
with f2 = 0. With the initial state (20), the time-dependent wave function is given by the 
expansion 



where £ = h/(2I), and the coefficients ( m are given in Eq. (21). As displayed in Fig. 2, 
the distribution |Cm| 2 typically separates in a number of peaks centered at fn — 2nl, where 
n = 0, ±1, ±2, .. .which are separated by 21. Thus, the superposition state (22) can be 



V. FREE EVOLUTION ON A RING 




m 



(22) 
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considered as a series of elementary wave packets centered at 2nl, in the angular-momentum 
space. Each of these peaks gives a separate contribution to the wave function that moves 
with its own angular group velocity 2£m = 4£n/ = v n . The angular separation between 
neighboring wavepackets is given by 4£lt, which is proportional to I. Since wave packets 
with opposite angular-momentum values will move in opposite directions, i. e. clockwise 
and anticlockwise, they will eventually meet again at some time t = T and start to interfere. 

In order to estimate the time value that interference sets in, we use the fact that for 
not too small arguments VLt/2 the Bessel function J„(f2r/2) with the maximal value is 
the one with n = n max — Vtr/2. Hence, the meeting time of the pair of the strongest 
counterpropagating packets is 



7T 



7T 



v max 2£ttrl 

The exact expression for the time-dependent wave function can be expressed in an integral 
form by using the mathematical identity [9] 

1 



exp (imcf) — i£tm 2 ) 



d<(> e im<p exp 



(0-0') /4# 



(23) 



which can be checked by performing the integration. When substituting this identity in the 
right-hand side of Eq. (22), and using the expansion (20), we arrive at the exact expression 

1 



$(0,t) 



d<t> $(0',O)exp 



10- 



/4# 



(24) 



J c 

A similar equation is well-known to describe the free evolution of a quantum particle in one 
dimension. In the present case it is crucial that the integration be performed over all values 
of (/>', while using that the wave function $((/>', 0) is periodic. Because of this periodicity, we 
can express the integral in (24) as a sum of bounded integrals 

r.27r(p+l) 



d>(<M) = -== Yl / d( f> $ (^0)exp iU-<t>) /^t 



(25) 



By a shift of variables the integrations can be performed over the interval [0, 2tt], which leads 
to an integral expression over a single interval 



$(<M) 



J2 exp[2(0-27rp) 2 /4^] 



p=— oo 



2tt 



, t ) exp 



2irp) 4' /2£t . 

(26) 



Here we introduced the modified wave function <3> which is just the initial wave function, 
modified by a phase factor, defined by 

(27) 



$ (j,\t) = $(0', 0) exp i(j>' 2 /4£f 



In order to emphasize its physical significance, we write Eq. (26) in the form 



$(0,t) 



1 



exp[i(<f>-2irp) 2 /4tt]F 



2rcp 



2& 



(28) 



p=— oo 

where the function F is the Fourier transform of the modified wave function defined over a 



single period 



exp 



(29) 



For a freely evolving quantum particle in one dimension, the time-dependent wave func- 
tion has the same form as the term with p = in (28). The other terms can be understood 
from the periodic nature of the dynamics on the circle, where each period of the initial 
wave function serves as an additional source that contributes to the wave function ^ (0, t) 
in the relevant interval [0, 2n]. Because of the finite range of the integration in Eq. (29), 
the distinction between the modified wave function and the initial wave function vanishes 
for times t obeying the inequality t ^> l/(£iV), when we find in a good approximation 



,t 



>',o). 



(30) 



In this limit, the function F is just the Fourier transform of the initial wave function $(0, 0), 
and $(0, t) is simply determined by the Fourier transform F of the initial wave function 
$(0, 0) multiplied by a phase factor. The equation (28) has the flavor of the far-field picture 
of the time-dependent wave function. The Fourier transform of the initial state determines 
not only the momentum wave function, but also the asymptotic form of the coordinate wave 
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function, scaled by a factor that varies linearly with time. Characteristic for the present 
case of evolution on a circle is that each interval of length 2ir serves as a separate source, 
each giving a contribution to $(0, t). Since the Fourier transform of the wave function 
determines the angular-momentum amplitudes, we may conclude that the wave function for 
not too small times has the same form as the initial distribution of angular momentum, 
scaled by the factor 2£t. 

It is clarifying to follow the temporal evolution of |$ (<p,t)\ 2 by distinguishing two time 
regions, namely < t < T and t > T. In the region < t < T, the wave function has not 
yet spread beyond a single period of length 2n, and only a single term in eq. (28) (or (26)) 
differs from zero. The contribution to the wave function coming from different sources do not 
overlap yet, so that one can neglect the interference term between them. At later times t > T, 
the diffraction pattern on the interval [0, 27r] is formed as an interference pattern between 
two and more terms in the superposition state (25). This picture is confirmed by numerical 
calculation of the diffraction pattern for the two time regimes. In Fig. 3 the angular 
probability distribution |$ (4>,t)\ 2 is shown for a time t <T. The spatial pattern resembles 
the angular momentum distribution shown in Fig. 2. Figure 4 displays the same probability 
distribution for a later time t > T. One notices that the counterrotating components give 
rise to clear interference fringes. These fringes will be quite sensitive to any perturbation in 
one of the arms. This suggests to use the present scheme as an atomic interferometer [9]. 
Usually, interferometers have two key components, namely a beam splitter and a mirror. A 
coherent incoming atomic beam is split into spatially separated components by the beam 
splitter. Two arms are getting formed, which freely propagate and may undergo different 
phase shifts, which are probed by recombining the two arms. The interference pattern 
contains the information of the phase perturbation in one of the arms. Recombination 
usually requires atomic mirrors. In atom optics, beam splitters and mirrors are commonly 
realized by using light pulses, with carefully selected duration and shape. 

In the present case, only a single pulse is required that splits the initial atomic wave 
packet into components rotating to the left and to the right. No mirrors are employed in 
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this scheme. Instead, one uses the radial potential U (p), to constrain the atomic motion to 
a ring. Radial potentials can be realized by hollow light beams, which are widely used in 
atomic interferometric schemes. 

VI. RADIAL DISPERSION 

In this chapter we consider the radial dynamics of the diffracted wave function during its 
free evolution, after the passage of the circular lattice. We assume that the wave function 
at time t = 0, after the diffracting pulse, is factorized as 

M>(p,0,O) = Q(p,O)$(0,O), (31) 

where the radial part Q of the wave function is sharply peaked at p = p , and the 
angular wave function is specified by Eq. (20). The radial function Q is normalized 
(J °° dpQ 2 (p)p = 1). We wish to study the possible deformation of the wave packet, when 
the radial dispersion is included during the stage of free evolution. We take the simplest 
possible trapping potential, which allows radial dispersion, and we take for the confining 
potential U (p) an infinitely deep cylindrical box with radius a, as defined by 

U (p) = for p < a, and U (p) = oo for p > a. (32) 

This potential models a hollow light beam. With this potential, the normalized eigenfunc- 
tions of the Hamiltonian for the cylindrical coordinates during the free-evolution stage take 
the form 

^ nm (p, 0) = R nm (p)^—e m \ (33) 
where the radial functions R nm are solutions of the equation 

Rnm(p) = £nm,Rnm(p), (34) 

with S nm are the corresponding eigenenergies. The radial functions are proportional to the 
Bessel function of order m 
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K 1 I d d 



2M 



p dp^ 'dp 



m 
~7z 



+ U(p) 



Rnm(p) oc J m {a nm p/a), (35) 

with R nm normalized in the interval < p < a. In order that the wave function vanishes at 
the edge p = a of the cylindrical well, we have to take the numbers a nm for various values 
of n as the subsequent zero's of the Bessel function J m . This determines the corresponding 
eigenenergies as 

£ nrn = h\a 2 nm . (36) 

with A = h/(2Ma 2 ). For each value of the angular momentum m, the set of functions 
Rnm(p) is complete. An expansion of the initial state (31) in the energy eigenfunction is 
found when we expand the initial radial wave function Q(p, 0) in the radial eigenfunctions 
(35), so that 

Q(p, 0) = ^ C nmRnm(p), (37) 

n 

while substituting Eq. (20) for the initial angular state $(0,0). For the time- dependent 
state we find 

*(p,0,t) = Y,-7K=(rne im(P Qrn(P,t), (38) 



/2tt 

in 

where the m-dependent radial wave function Q m is 

Q m (p,t) = ^2c nm R nm (p) exp(-i£ nm t/h). (39) 

n 

/,From Eq. (37) one notices that Q m (p,t) = Q(p,t), independent of the angular momentum 
m. It is obvious from the radial Schrodinger equation (34) and the initial condition (31) 
that the normalized radial wave function obeys the identity Q m (p,t) = Q- m (p,t) for all 
to. Moreover, since the total wave function before diffraction is even in 0, it must remain 
even for all times. This implies that ( m = £_ m for all to. So just as discussed in Sec. IV, 
the angular distribution separates in different wave packets that are counterrotating. Since 
the phase of ( m Q m is even in to, its derivative with respect to to will be odd, and the 
angular group velocities of packets with opposite values of to will be opposite. This leads 
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to interference after the packets have traversed the entire ring. The initial radial function is 
taken as a narrow Gaussian 



Here L is the width and p represents the initial position of the wave packet within the 
box. The normalized wave function Q m (p,t) describes the radial dynamics for each value 
of the angular momentum m. As an example, we evaluate the time behavior of the average 
radius for each angular momentum, with the given initial radial state (40), according to the 
expression 



The result is displayed in Fig. 5, in the special case that m = 10. The average radius 
displays oscillations, which can be understood as arising from the outward motion due 
to the centrifugal potential, followed by reflection at the hard wall of the cylinder. The 
oscillations display collapse, followed by a revival. These may be viewed as arising from the 
initial dephasing of the contribution from the radial eigenfunctions R nm with different values 
of n, due to their energy difference. The revival of the oscillation can be understood from 
the discrete nature of the contributing energy eigenvalues, when the phase factors due to 
neighboring eigenenergies have built up a phase difference 2n. Because of the conservation of 
angular momentum, the probability density near the origin remains zero. The interference 
between the counterrotating wave packets is illustrated in Fig. 6, for the ring at radius 
P — Po — a /2- Fig. 6a shows the short-time separation of the angular wave packets. Fig. 
6b displays the interference that arises as soon as overlap occurs around = re between 
the clockwise and the anti-clockwise rotating packets. This demonstrates that the radial 
wave functions Q m (p,t) for different values of m have sufficient overlap, so that the angular 
interference survives the effect of radial dispersion. 




(40) 




rn 
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VII. CONCLUSIONS 



In this paper we describe the diffraction of an atomic wave by a circular optical lattice. 
Such a lattice can be formed by the superposition of two Laguerre-Gaussian beams with 
opposite helicity, which gives rise to a standing wave in the angular direction. Such a light 
field will split a single localized atom into clockwise and anticlockwise rotating components. 
If the system is in a trapping potential in the form of a ring or in a cylindrical box, these 
counterrotating components give rise to interference. We express the spatial pattern in the 
far diffraction field in terms of the Fourier transform of the near-field diffraction pattern. 
The periodic nature of the circular motion modifies this relation compared with the case of 
diffraction by a linear standing wave. The general conclusions are backed up by numerical 
calculations. Characteristic for the circular case is that the wave packets corresponding to 
opposite angular momentum will cross each other, even without applying light pulses to 
reverse their motion, as in more common interferometric schemes. The scheme is reasonably 
robust to changes in the radial confining potential. 
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X 

FIG. 1. Circular lattice structure due to the trapping potential V(x,y). The plot shows 
V (x, y) /Ml for I = 2, 4 in Cartesian coordinates. 
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-100 -50 50 100 

FIG. 2. Probability distribution of angular momentum m before (upper) and after the pulse 

(lower). Here the helicity of the circular lattice is I = 10, the initial state is determined by N = 10 

and the pulse duration t is given by S7r = 6. 
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-71 -Till 7T/2 71 

FIG. 3. Angular distribution | i)| 2 is plotted versus the azimuthal angle (f> before the left 
and right rotating components cross. Here £t = 3tt x 10~ 3 , the value of N determining the width 
of the initial state, the helicity I and the pulse duration r are the same as in Fig. 2. 
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FIG. 4. Angular distribution | t)\ 2 is plotted versus the azimuthal angle <fi after the left and 

right rotating components cross. Here £i = 6ir x 10 -3 , the value of N, the helicity / and the pulse 

duration r are the same as in Fig. 2. 
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0.1 0.2 0.3 0.4 ^ t 

FIG. 5. Time behavior of the average radial distance {pit)) /a for the angular momentum 
m = 10. N is the same as in Fig. 2, the width of the initial Gaussian is L = 0.01a, and the initial 
average radial distance is po = a/2. 
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7t/2 n 37i/2 ± 2% 

FIG. 6. Angular distribution \^ (po, (f),t)\ 2 in the presence of radial dispersion, at the ring 
P — Po = o/2- The time values are determined by Aii = ir x 10 -3 and Xt2 = 2tt x 10~ 3 . N,1,Qt 
are the same as in Fig. 2. 
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